j211filtfilt 
% j211Filtfilt.m 
%*4 Description 
j211filtfilt -- Zero-phase filtering method that conforms to SAE J211 


SEE: SAE J211-1 for more information 


ALSO SEE <j211_response.html |j211_response|> for frequency and phase plots of the 
ilters 

9 
Note that for the sampling rate you are using, you will need to make sure 
that the filter meets the corridors specified in J211. We have found that 
for sampling rates below ~15,@@@ Hz, the original scale factor given 
(2.0775) is insufficient to meet the corridors. use j211 response.m to 
check filter response before use. This code is provided without warranty. 
Use at your own risk. 
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% Usage 
% y = j211filtfilt(cfc,samplerate,x) 
%*4 Input Arguments 


% * *cfc* | |double| | _cfc_ is the filter class (can be 6@, 108, 180, 600, or 1000 
Hz) 

% * *samplerate* | |double| | _samplerate_ is the sampling rate of the 

% input data array (_x_) 

% * *x* | |array| | _x_ is the input data array to be filtered 

4% Output Arguments 

% * *y* | |array| | _y_ is the filtered data 
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Revision Log: 

Version Date Modified By Reason 

2.0 21-MAY-@9 J. Somers Wrote to work with my scripts and to 
make sure the filter falls within 
the SAE corridor. Added padding to 
the beginning and end to minimize 
startup effects 

2.1 22-SEP-14 J. Somers Added b and a outputs 

2.1.1 @3-JUN-17 J. Somers Formatted for commercial crew 
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% Code 
function [y,b,a] = j211filtfilt(cfc,samplerate, x) 
hh 
% *Check Input Arguments* 
if size(x,1) > size(x,2) 
XS 5C- Ss 
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flag 
else 

flag 
end 


i] 
= 
Noe 


8; 


len = length(x) ; 


hh 


% *Setup filter specs* 
T=1/samplerate; %sample period in seconds 
% Sf = 2.0775 ; % actual scale factor from J211. Smpf dependent 


switch cfc 
case 


requirements 


case 


requirements 


case 


requirements 


case 


requirements 


case 


requirements 


60, 


sf = 2.2 ; 


d 
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if samplerate < 2700 
error([ ‘Filter 


end 
108, 
sf = 2.2 


3 


"for a 


if samplerate < 5000 
error([ ‘Filter 


end 
180, 


sf = 2.2 ; 


d 


"for a 


if samplerate < 6200 
error([ ‘Filter 


end 
600, 


sf = 2.0 ; 


J 


"for a 


if samplerate < 7000 
error([ ‘Filter 


end 
100ee, 


"for a 


if samplerate < 9600 
error([ ‘Filter 


end 


"for a 


if samplerate < 6442*2 


else 


sf = 2.3 5 


d 


Performance does not meet SAE 


num2str(samplerate) ' Hz']) 


Performance does not meet SAE 


num2str(samplerate) ' Hz']) 


Performance does not meet SAE 


num2str(samplerate) ' Hz']) 


Performance does not meet SAE 


num2str(samplerate) ' Hz']) 


Performance does not meet SAE 


num2str(samplerate) ' Hz']) 
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J211 


J211 


J211 


J211 


J211 
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Sfos 2-3 
end 
otherwise, 
sf = 2.2 ; 


end 

wd=2*pitcfc*sf ; %cfc = channel filter class 
wa=sin(wd*T/2)/cos(wd*T/2); 
aQ@=wa*2/(1.0+sqrt(2)*watwa’2) ; 

al=2*aQ; 

a2=aQ; 

be=1; 

b1=-2* (wa*2-1)/(1+sqrt(2)*watwa’2) ; 
b2=(-1+sqrt(2)*wa-wa*2) /(1+sqrt(2) *watwa%2) ; 

a=[b@ -b1 -b2]; %coefficients for filter from SAE J211 
b=[a@ al a2]; %coefficients for filter from SAE J211 


wh 

% *Pad data at beginning and end* 

numPts = floor(1@ * samplerate / 1000) ; %% 10 ms of data 
x1 = zeros(1,numPts); 

x = [x1 x x1] ; 


%% 
% *Filter data* 
y=filtfilt(b,a,x); 


hh 

% *Remove padding* 

if length(numPts:length(y)-numPts) > len 
y = y(numPts+1:end-numPts) ; 


else 
y = y(numPts:end-numPts) ; 
end 
if flag 
yey 
end 
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